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1 Introduction 



The purpose of Pegase is the study of galaxies by evolutionary synthesis. This version supersedes 
our previous model (Fioc & Rocca-Volmerange 1997; contributions in Leitherer et al. 1996). The main 
differences are the implementation of: 

• stellar evolutionary tracks with non-solar metallicities; 

• the library of stellar spectra of Lejeune et al. (1997, 1998); 

• radiative transfer computations to model the extinction. 

The extension to the far- infrared (Fioc & Dwek, in prep.) and a detailed modeling of the nebular emission 
(Moy, Rocca-Volmerange & Fioc, in prep.) are in progress. Synthetic spectra computed from standard 
star formation scenarios fitted on new statistical templates for nearby galaxies (Fioc & Rocca-Volmerange 
1999; Fioc & Rocca-Volmerange, in prep.) will be proposed in a near future, as well as their colors, k- 
and e-corrections. 

To be informed of the future developments of Pegase, to require specific computations, ask questions, 
or make comments or suggestions, mail us at pegase@iap.fr. 

2 Contents of the directory 

2.1 List of files 

README.tex 
README . ps 
SSPs.f 
calib . f 
colors . f 
scenarios . f 
spectra. f 
IMF_Kennicutt.dat 
IMF_Kroupa.dat 
IMF_MillerScalo . dat 
IMF_Salpeter.dat 
IMF_Scalo86.dat 

2.2 Codes 

• calib.f : code computing the calibrations of the filters. 

• colors. f: code computing colors and other quantities. 

• spectra. f: code computing synthetic spectra of galaxies and other quantities. 

• SSPs.f: code computing the properties of simple stellar populations (SSPs), i.e. populations of 
stars formed simultaneously with the same metallicity. 

• scenarios. f: code used to prepare the input file (star formation scenarios) to spectra. 



IMF_Scalo98.dat 
Spitzer . dat 
WW . dat 
ages . dat 
calib . dat 
dust . dat 
filters . dat 
list_IMFs.dat 
list_tracks . dat 
HI I. dat 

BD+17o4708.dat 
stellibCM.dat 



stellibLCBcor . dat 
SunLCB . dat 
VegaLCB . dat 
King. dat 
slab . dat 

tracksZ0.0001.dat 
tracksZO . 0004 . dat 
tracksZ0.004.dat 
tracksZO . 008 . dat 
tracksZ0.02.dat 
tracksZO . 05 . dat 
tracksZO. 1. dat 
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The codes are written in Fortran 77. Though available on most systems, some features are non-standard: 

• use of lowercase letters and underscore; 

• use of identifiers longer than 6 characters; 

• implicit none; 

• do while; 

• do ... end do; 

• list-directed input /output in internal files. 
Fortran 90 should also work. 

To compile and execute the file name.f (name=calib/colors/spectra/SSPs/scenarios): 

Unix: 

• f 77 name . f -o name 

• name 

VMS: 

$ for name 
$ link name 
$ run name 

You may have to rename name.f as name. for and change the lowercase letters to uppercase. 
2.3 Data files 

2.3.1 Stellar evolutionary tracks 

Stellar evolutionary tracks for various metallicities (Z) and helium abundances (Y): 

• tracksZ0.0001.dat: Z = 0.0001, Y = 0.23. 

• tracksZ0.0004.dat: Z = 0.0004, Y = 0.23. 

• tracksZ0.004.dat: Z = 0.004, Y = 0.24. 

• tracksZ0.008.dat: Z = 0.008, Y = 0.25. 

• tracksZ0.02.dat: Z = 0.02, Y = 0.28. 

• tracksZ0.05.dat: Z = 0.05, Y = 0.35. 

• tracksZO.l.dat: Z = 0.1, Y = 0.48. 
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The names of these files are written in listjtracks . dat (read by SSPs). 

The tracks proposed here come mainly from the "Padova" group. At Z = 0.1, pseudo-tracks for masses 
larger than 9 Mq have been computed from the corresponding masses in the Z = 0.02 and Z = 0.05 
sets. For stars undergoing the helium flash, the zero-age main sequence tracks are connected to the 
zero-age horizontal branch tracks assuming a Reimers law for the mass loss along the first giant branch 
with rj = 0.4. The same law is used to describe the mass loss during the early asymptotic giant branch 
phase (EAGB). Pseudo-tracks are then computed for the thermally pulsing AGB (TPAGB) phase using 
the equations proposed by Groenewegen & de Jong (1993) with rj = 4 (van den Hoek & Groenewegen 
1997). Hydrogen burning post- AGB and CO white dwarf tracks from Blocker (1995) and Schonberner 
(1983) are then connected. For low-mass stars becoming helium white dwarfs, but for which the Padova 
group does not provide the tracks, we use the Althaus & Benvenuto (1997) models. The positions of 
(unevolving) low-mass stars in the HR diagram come from Chabrier & Baraffe (1997). 

2.3.2 Initial mass functions 

Some initial mass functions (IMF) are already defined analytically in SSPs.f: 

• In: lognormal IMF (Miller & Scalo 1979). 

• RB: Rana k Basu (1992). 

• Fe: Ferrini (1991). 

Others are given in specific files: 

• IMF_Kennicutt.dat: Kennicutt (1983). 

• IMF_Kroupa.dat: Kroupa et al. (1993). 

• IMF_MillerScalo.dat: Miller & Scalo (1979). 

• IMF_Salpeter.dat: Salpeter (1955). 

• IMF_Scalo86.dat: Scalo (1986). 

• IMF_Scalo98.dat: Scalo (1998). 

These files are listed in list_IMFs.dat. 

2.3.3 Filters and calibrations 

Filters: The passbands of the filters are provided in the file filters.dat. Since you may want to 
change it to add new filters, we detail its content here. The structure of filters.dat is the following: 

• 1 st line: number of filters (Agiters)- 

• -^filters blocks, one for each filter, containing: 

- 1 st line: 

* the number of wavelengths (AT wave i engths ); 

* the type of transmission (see below); 
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* the type of calibration (see below); 

* a code between quotes ('...') identifying the filter and used in calib.dat; 

* reference, comments (optional). 
- ^"wavelengths hnes containing: 

* the wavelength in A; 

* the transmission at this wavelength. 

Type of transmission: 

• 0: the shape of the transmission curve (T\ = T u ) corresponds to the energy transmitted. 

• 1: the shape of the transmission curve corresponds to the number of photons transmitted. 

• 2: used for 

r 4250A 

F v dA 

3950A 



Atooo = (Bruzual 1983). 



F v dA 

3750A 



Type of calibration: 

• 0: used for -D4ooo- 

• 1: standard system 



f F\(*) T\d\ 

m(*) = — 2.51og 10 —d h 0.03 (i.e., the magnitude of Vega is 0.03). 

' J F \ (Vega) T\ dA 

• 2: AB system; used for the SDSS filters (u', g', r', i', z') 

f F u (*)T u &v 

m A B(*) = -2.51og 10 48.60 {F u in erg.s _1 .cm _2 .Hz _1 ). 

J T u dv 

• 3: Thuan & Gunn system; used for the Thuan & Gunn filters (u, v, g, r) 

f F\ (*) T\ dA 

titgW = -2.51og 10 — r + 9.50. 

/ F A (BD + 17°4708)T A dA 

• 4: used for the WFPC2 filters (F300W, F450W, FW6W, F8UW) and the ultraviolet filters at 
1650 A, 2500 A and 3150 A 

f F\ (*) T\ dA 

"i-2i 10 (*) = -2.51og 10 ± j. 21.10 (Fx in erg.s _1 .cm" 2 .A _1 ). 

J T x d\ 
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• 5: used for the FOCA filter at 2000 A 



"1-21.175 (*) = -2.51og 10 



/ 



F \ (*) T\ dA 



-21.175 



(F\ in erg 



.s _1 .cm _2 .A _1 ). 




Calibrations: The calibrations of the filters are computed by the code calib.f and written in 
calib.dat. The structure of this file is the following: 

• 1 st line: caption of the file. 

• One line for each filter containing: 

— the name of the filter; 

— the corresponding index used in colors; 

— the apparent flux of Vega in erg.s - 1 .cm — 2 .A — 1 : / Fa (Vega) T\ dA/ / T\ dA; 

— the "area" of the filter in A: J T\ dA; 

— the mean wavelength in A: A = J A T\ dA/ J T\ dA; 

— the effective wavelength of Vega in A: J X Fx (Vega) Tx dA/ / Fx (Vega) T\ dA; 

— the AB-magnitude of Vega; 

— the Thuan & Gunn-magnitude of Vega (99.999 if undefined); 

— the "monochromatic" luminosity of the Sun in erg.s -1 . A -1 : J F\(Q) T\ dA/ jT^dA. 

The AB-magnitude toab may be computed from the standard magnitude m (in the Vega-system) as: 
"^abM = m(*) + mAB(Vega) and the same for the Thuan & Gunn-magnitude. You may also directly 
modify the type of calibration in filters.dat to change the default ones used in colors. f. 

2.3.4 Other files 

• stellibCMcor.dat: stellar library of Clegg & Middlemass (1987); T cS > 50000 K. 

• stellibLCBcor.dat: stellar library of Lejeune et al. (1997, 1998; corrected version (BaSeL-2.0)); 
T cS < 50000 K.. 

• ages .dat: ages at which the synthetic spectra will be written. 

• HI I. dat: used to compute the nebular emission (continuum and lines). 

• dust . dat: extinction properties of graphites and silicates (Draine & Lee 1993; Laor & Draine 1993). 

• slab. dat: results of the radiative transfer code for an homogeneous slab model for both stars and 
dust (Fioc 1997); used to model the extinction for disk galaxies. 

• King. dat: results of the radiative transfer code for a spheroidal geometry, where the stars are 
distributed according to a King profile and the dust to a power | of the King profile (Fioc Sz 
Rocca-Volmerange 1997); used to model the extinction for elliptical galaxies. 

• VegaLCB.dat: spectrum of Vega (Lejeune et al. 1997; computed by R.L. Kurucz). 
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• BD+17o4708.dat: spectrum of the F subdwarf BD+17°4708 (Oke & Gunn 1983) used to calibrate 
the Thuan & Gunn (1976) photometric system. 

• SunLCB.dat: spectrum of the Sun (Lejeune et al. 1997; computed by R.L. Kurucz). 

• Spitzer.dat: table 5.4 of Spitzer (1978, p. 113). 

• WW.dat: stellar yields of Woosley & Weather (1995). 

3 Computing synthetic spectra 
3.1 Preliminaries 

Stars [> Except for the star formation rate, which takes also into account substellar objects, "star", 
"stellar", etc., refer only to luminous stars to the exclusion of stellar remnants (old white dwarfs, 
neutrons stars and black holes) and substellar objects. 

Gas > Gas means both the gas strictly speaking and the dust. 

Metallicity > All the metallicities are given in mass fraction. 

Galaxy, reservoir... [> We consider only the baryonic matter (with the constant mass M tot ) and dis- 
tinguish two zones: 

• The galaxy itself (mass M ga i). Unless otherwise specified, all the quantities refer only to this 
zone. 

• A reservoir of gas only surrounding the galaxy (mass M res ). 
Initially, both zones contain only gas and we have either 

• M g3i \ = M to t and M rcs = 0: the galaxy is already fully constituted; 

or 

• Mg a i = and M res = M to t'- the galaxy forms entirely by infall from the reservoir. 

In both cases, the reservoir may be replenished by galactic winds occurring in the galaxy. These 
moreover interrupt the infall. 

Normalized quantities [f, J] > Some quantities in the following are normalized to M to t = 1M@. To 
obtain the value for a given M to t, you have either to: 

• multiply them by M tot [in M Q ]: quantities denoted by a "f"; 

or 

• to add — 2.5 log 10 M tot [in Mq]: quantities denoted by a "J". 

Maximal star formation rate [> If, at any time t, the normalized star formation rate SFR(t) exceeds 
SFR mSLX , its maximal possible value given the amount of gas available, spectra sets SFR(t) to 
SFR max and, the first time it happens, also prints a warning on the screen and in the header of the 
output file. 
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Inclination-dependent quantities [§] > If there is some extinction in the disk geometry, the emission 
is not isotropic. If you choose to compute the spectra for a specific inclination, the monochromatic 
luminosity L\ is then defined as L\{Qq) = 4irA\(6 ), where A^(6*o) dAdw(^o) is the energy radiated 
between A and A + dA and escaping from the galaxy in a solid angle du(6) = 2irsm9d9 having an 
inclination 9 = O to the axis of rotational symmetry. Inclination-dependent quantities (monochro- 
matic or in-line luminosities, magnitudes, etc.) are denoted by a "§" in the following. This does 
not apply to the bolometric luminosity, as computed here (Lbol = jj Aa(0) d\did(6) ^ J L\(9q) dA), 
nor to the dust emission, which is supposed to be isotropic (negligible self-absorption). You may 
also output inclination-averaged spectra (La = / Aa(0) did (9)) rather than for a specific inclination. 

Consistent evolution of the metallicity > Stars are formed with the same metallicity as the ISM. 
3.2 Procedure 

The synthetic spectra are computed in three steps: 

1. run SSPs 1 to compute the properties of SSPs of different metallicities; 

2. run scenarios to prepare the input file to spectra containing the parameters of the star formation 
scenarios; 

3. run spectra. 



3.2.1 SSPs 

You will be asked: 

• The shape of the initial mass function (enter the corresponding number). 

• The lower mass of the IMF. 

• The upper mass. 

• The type of supernovae ejecta. 

• If you want to take into account the ejecta due to stellar winds in high-mass stars through a 
somewhat dubious procedure. 

• A prefix (e.g. prefix). The output files corresponding to the tracks tracksZ* . dat will be named 
pre/ia:_tracksZ*.dat 2 and will be listed in the file called prefix_SSPs . dat. 

Default values > Default values are proposed for some quantities. Type <return> to select them. 

You do not need to run SSPs every time if you keep the same IMF and the other parameters asked by SSPs. 
2 spectra will interpolate between the resulting files. If the metallicity of the stars formed in spectra is lower (resp. 
higher) than the lowest (resp. highest) metallicity of every file, spectra does not extrapolate but uses the data of the file 
with the lowest (resp. highest) metallicity. 
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3.2.2 scenarios 
You will be asked: 

• the name of the output file, e.g. scenarios . dat (it must be a new name); 

• the name of the file (prefixSSPs . dat in the example above) listing the names of the SSP files; 

• the fraction of close binary systems (this quantity is used to compute the number and the ejecta of 
SNIa, assuming the W7 model of Thielemann et al. (1986) and the formalism of Greggio & Renzini 
(1983) and Matteucci & Greggio (1986)). 

These data are common to all the star formation scenarios chosen later. 
Then for each scenario, you will be asked: 

• The name of the file containing the corresponding synthetic spectra (just type end to stop). 

• The initial metallicity of the interstellar medium (ISM). 

• Whether you want to build your galaxy by infall or prefer to start from a galaxy already constituted. 
The infall rate, computed as a function of the time t, is: 

,j exp(-t/t infa ii) 

■Mtot 7 • 

Hnfall 

You will have to provide ii n f a ii (Myr) and the metallicity of the infalling gas. 

• The type of star formation scenario (characterized by an integer) giving SFR [f] - the normalized 
star formation rate in M .Myr _1 - as a function of the time in Myr, the normalized mass of gas 
M gas [f] in M© and other quantities. 

- Types to 9 are reserved for predefined laws of star formation implemented in spectra. f: 

* 0: instantaneous burst: SFR(t) = 5(t). 

* 1: constant star formation rate: 

SFR{t)= Pl ift<j>2, 
= if t > p 2 . 
[pi] =M .Myr- 1 ; [p 2 ] = Myr. 

* 2: exponentially decreasing or increasing star formation rate: 

SFR(t) = p 2 eM - t/Pl) . 

Pi 

[ Pl ] = Myr; [p 2 ] = M . 

* 3: star formation rate proportional to some power of the mass of gas: 



SFR(t) = 



P2 

[pi] = 1; N =Myr.M" 1 . 
* 4 ... 9: not yet defined. 
You will then be asked the values of the parameters (pi, p 2 , etc.). They must be real. 

Types >10: you have to implement your star formation law in spectra. f (see section 5.2). 
You will be asked the number of parameters used by this law and the values (real) of each 
one. 
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— Types -1 and -2 are for files containing the star formation rate as a function of time. You 
will be asked the name of the file (e.g. SFRfile): 

* -1: SFRfile must contain on each line the age in Myr and SFR separated by blanks. 

* -2: SFRfile must contain on each line the age in Myr, SFR and the metallicity of the 
forming stars separated by blanks. This metallicity may be inconsistent with that of the 
ISM. 

These quantities must be real. The first age in SFRfile must be 0. and the last must be higher 
than 20000. The computation of the star formation rate at intermediate ages is performed by 
spectra. 

• If the type of star formation scenario is not -2, whether you want a consistent evolution or prefer 
to form stars with a constant metallicity (asked later). 

• The fraction (in mass) of the star formation rate used to form substellar objects. These objects 
lock the mass and are supposed to emit no light. 

• If you want galactic winds. Galactic winds expel all the interstellar medium from the galaxy after 
a given time (asked later) and prevent any further star formation. 

• If you want to take into account the nebular emission, i.e. the continuum and lines emitted by 
the ionized gas in star- forming regions. The emission in the continuum and the hydrogen lines is 
computed from the number of Lyman continuum photons in the case B of recombination. Typical 
observed ratios to H/? are taken for other lines. 

If you hereafter choose to have some extinction, a fraction of the Lyman continuum photons will 
be absorbed by the dust inside the Hn region rather than by the gas. This fraction is computed 
according to the prescriptions of Spitzer (1978, p. 113) and assuming that 70% of the Lyman 
continuum photons are absorbed by the gas at solar metallicity. 

• If you want to introduce some extinction: 

— 0: No extinction. 

— 1: Extinction for a spheroidal geometry. 

— 2: Extinction for a disk geometry; inclination-averaged. 

— 3: Extinction for a disk geometry; specific inclination. You will then be asked the inclination 
in degrees relative to face-on. 

The optical depth is estimated from the mass of gas and the metallicity. The absorption, the 
albedo and the asymmetry parameter are computed from Draine & Lee (1984) and Laor & Draine 
(1993) data for a mixture of graphites and silicates depending on the metallicity and fitted on the 
Magellanic Clouds and the Milky Way (cf. Pei (1992)). 

In the cases 1 and 2, all the Lyman continuum photons not absorbed by the gas as well as those 
emitted in the Lya line are absorbed by the dust as soon as the metallicity of the ISM is non 0. 

Default values > Default answers are proposed for some questions. Just type <return> to select them. 
Default names of the output files are created by inserting the number of the scenario between 
the prefix spectra and the suffix .dat (see however note 3). For the other questions, the default 
answers are those defined in scenarios. f the first time you answer to a specific question. When 
you have already answered to this question for a previous scenario, the default is your last choice. 
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3.2.3 spectra 

Type the name of the file of scenarios (scenarios . dat in the example above) when required 3 . 
The structure of the output files is the following: 

• A block describing the evolutionary scenario and ending with a line of asterisks (*** . . . ***) 
only. 

• One line with: 

— the number of timesteps (A^i mesteps ); 

— the number of wavelengths of the continuum (A con tmuum); 

— the number of emission lines (N\\ nes ). 

• A block containing the -/V cont i nuum wavelengths (A) of the continuum (5 per line). 

• A block containing the N\ [ncs wavelengths (A) of the emission lines (5 per line). 

• Ati mest eps blocks (one for each timestep) containing: 

— 1 st line: 

* the time (Myr, integer); 

* the normalized mass of the galaxy [f] (Mq); 

* the normalized mass in stars [f] (Mq); 

* the normalized mass in white dwarfs [f] (Mq); 

* the normalized mass in neutron stars and black holes [f] (Mq); 

* the normalized mass in substellar objects [f] (Mq); 

* the normalized mass in the gas [f] (Mq); 

* the metallicity of the interstellar medium (mass fraction) ; 

* the mean metallicity of stars averaged on the mass (i.e., the mean initial metallicity of the 
stars still alive averaged on their initial mass); 

* the mean metallicity of stars averaged on the bolometric luminosity (i.e., the mean initial 
metallicity of the stars still alive averaged on their present bolometric luminosity). 

— 2 nd line: 

* the normalized bolometric luminosity [f] (erg.s -1 ); 

* the optical depth in the V-band (5500 A) from side to side (through the center for the 
spheroidal geometry and along the axis of rotational symmetry for the disk geometry); 

* the ratio of the luminosity emitted by the dust to the bolometric luminosity; 

* the normalized star formation rate [f] (M Q .Myr _1 ); 

* the normalized number of Lyman continuum photons emitted [|] (s -1 ); 

* the normalized SNII rate [|] (Myr -1 ); 

* the normalized SNIa rate [f] (Myr -1 ); 

* the mean age of the stars averaged on the mass (Myr); 

* the mean age of stars averaged on the bolometric luminosity (Myr). 

3 If one of the files of spectra you want to create already exists, spectra appends one or more "+" to the name of the new 
file and prints a warning on the screen. 
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— A block containing the normalized monochromatic luminosities [f, §] (erg.s 1 .A 1 ) of the 
^continuum wavelengths of the continuum (5 per line). 

— A block containing the normalized luminosities [f, §] (erg.s -1 ) of the Ny mcs emission lines (5 
per line). 

4 Computing colors 

To compute colors, luminosities, etc., for a given set of spectra, run colors and type the name of the 
input file (spectra) when required. You are then asked the name of the output file (colors). If you just 
type <return>, the name of the output file is created by adding the prefix colors, to the name of the 
input file. 

The structure of the output file is the following: 

• A block describing the evolutionary scenario and ending with a line of asterisks (*** . . . ***) 
only. 

• One line giving the number of timesteps (iVtimesteps)- 

• Eight blocks consisting each in: 

— one line describing the quantity in each column; 

~ Atimcstcps lines giving these quantities 4 at each timestep. 

The quantities printed in the output file are the following: 

• 1 st block: time Mgal M* MWD MBHNS Mgas Zgas <Z*>mass <Z*>Lbol 

— time: time (Myr, integer). 

— Mgal [f]: normalized mass of the galaxy (M Q ). 

— M* [f]: normalized mass in stars (M Q ). 

— MWD [f]: normalized mass in white dwarfs (M ). 

— MBHNS [f]: normalized mass in neutron stars and black holes (M Q ). 

— Msub [f]: normalized mass in substellar objects (M Q ). 

— Mgas [f]: normalized mass in the gas (M Q ). 

— Zgas: metallicity of the gas. 

— <Z*>mass: mean stellar metallicity averaged on the mass. 

— <Z*>Lbol: mean stellar metallicity averaged on the bolometric luminosity. 

• 2 nd block: time Lbol tauV Ldust/Lbol SFR nSNII nSNIa <t*>mass <t*>Lbol 

— Lbol [f]: normalized bolometric luminosity (erg.s -1 ). 

— tauV: optical depth in the F-band. 

— Ldust/Lbol: ratio of the luminosity of the dust to the bolometric luminosity. 

— SFR [f]: normalized star formation rate (M .Myr -1 ). 

4 If no stars have formed yet, all the quantities are set to 0. This happens in particular at t = when the galaxy forms 
by infall. 
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— nSNII [f]: normalized rate of type II supernovae (Myr -1 ). 

— nSNIa [f]: normalized rate of type la supernovae (Myr -1 ). 

— <t*>mass: mean stellar age averaged on the mass (Myr). 

— <t*>Lbol: mean stellar age averaged on the bolometric luminosity (Myr). 

• 3 rd block: time nLymcont L(Ha) W(Ha) L(Hb) W(Hb) LB/LBsol LV/LVsol D4000 

— nLymcont [f]: normalized number of Lyman continuum photons emitted (s -1 ). 

— L(Ha) [f, §]: normalized luminosity of the emission line Ha (erg.s" 1 ). 

— W(Ha) [§]: equivalent width of the emission line Ha (A). 

— L(Hb) [f, §]: normalized luminosity of the emission line H/3 (erg.s -1 ). 

— W(Hb) [§]: equivalent width of the emission line H/3 (A). 

— LB/LBsol [f , §]: normalized blue 5 luminosity (Lb = J B L\ T\ dX/ J B T\ dA) in units of the solar 
blue luminosity, (i.e., L-q/L-q (©), which is different of X^L-q/Lq where L Q is the bolometric 
luminosity of the Sun). 

— LV/LVsol [f, §]: normalized visual luminosity (Ly/Ly(Q)). 

— D4000 [§]: intensity of the Balmer break (£>4ooo)- 

• 4 th block: time Mbol V U-B B-V V-K V-RC V-IC J-H H-K 

— Mbol [X]: normalized bolometric magnitude (Mb i(©) = 4.75). 

— V, U, B [$, §]: normalized absolute magnitudes in the filters of Buser & Kurucz (1978) [see 
note 5]. 

— RC and IC [X, §]: normalized absolute magnitudes in the R and / Cousins filters (Bessel 1990). 

— J, H and K, as well as L and M (see below) [X, §]: normalized absolute magnitudes in the filters 
of Bessel & Brett (1988). 

• 5 th block: time K-L L-M V-RJ V-IJ JK-V UK-JK JK-FK FK-NK 2000-V 

— RJ and I J [X, §]: normalized absolute magnitudes in the R and / Johnson filters (Johnson 
1965). 

— UK, NK [X, §]: normalized absolute magnitudes in the U and N filters of Koo (1986) 

— JK, FK [X, §]: normalized absolute magnitudes in the J and F filters of Kron (1980). 

— 2000 [X, §]: normalized absolute magnitude in the ultraviolet filter (Jose Donas, private com- 
munication) of the FOCA experiment (Milliard et al. 1991). 

• 6 th block: time V-ID ID-JD JD-KD BJ-V BJ-RF V-606 300-450 450-606 606-814 

— ID, JD, KD [X, §]: normalized absolute magnitudes in the /, J and K DENIS filters (Eric 
Copet, private communication); the passband of the K filter is the one determined at ambient 
temperature. 

— BJ, RF [X, §]: normalized absolute magnitudes in the Bj and Rp photographic filters (Couch 
& Newell 1980). 

For the sake of the consistency with the stellar library of Lejeune et al. (1997, 1998), the U, B, V filters used in the 
files of colors are from Buser & Kurucz (1978), not from Bessel (1990). The B filter is always B3, except for U — B where 
we use B2, which, as U3, is not corrected for the atmospheric absorption. 
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- 300, 450, 606, 814 [|, §]: normalized absolute magnitudes in the F300W, F450W, F606W, 
F814W filters of the WFPC2 instrument on the Hubble Space Telescope. 

• 7 th block: time u'-g' g'-r' V-r J r'-i' i'-z' u-v v-g g-V g-r 

- u' , g ' ,r ' , i ' , z ' [J, §] : normalized absolute magnitudes in the Sloan Digital Sky Survey filters 
(Fukugita et al. 1996). 

- u, v, g, r [$, §]: normalized absolute magnitudes in the Thuan & Gunn (1976) filters. 

• 8 th block: time 1650-B 1650-2500 3150-B 

1650, 2500 and 3150 [\, §]: normalized absolute magnitudes in Gaussian filters centered on 
the corresponding wavelengths in A of the Rifatto et al. (1995) data. 

Note that not all the filters provided in filters . dat are used in the ouput file of colors. See section 5.5 
if you want to use them. 

5 Adaptations 

5.1 IMF 

You may define your IMF as a series of p continuous piecewise power laws giving the number of stars n 
as a function of their mass m: 

dn 

if m£ ffli,mj + i , — oc m 1 (l<z<p). 

dlnm 

Create a file like this (see for example IMF_Scalo86.dat): 

p 

ra\ s\ 
m 2 s 2 

irip Sp 
m p+1 

and add its name at the end of the file list_IMFs . dat (type <return> at the end of the file). The lower 

mass should preferably be larger than 0.09 M Q and the upper mass less than 120 M Q to be in agreement 

rm p+ i d n 

with the tracks. The continuity of the power laws and the normalization of / - — m dm to 1 M & are 

J mi dm 

ensured by SSPs. 

5.2 Star formation rate 

You may define your own star formation rate. Search for the lines 
c if (typeSFR.eq.n>=10) then 

c SFR(i)=your SFR law (note that i = time in Myr + 1) 

c end if 

in spectra, f. Uncomment and modify them; then, express the normalized star formation rate [f] SFR(i) 
at the timestep i as a function of the age time (i) =i-l . , the normalized mass of gas [f] (sigmagas (i) ) or 
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other quantities. SFR(i) may also depend on free parameters param(l), param(2), ... , param(nparam) 
that you will have to provide when running scenarios. The maximal number of parameters nmaxparam 
is set to 99 in the declarations at the top of spectra. f; it should be enough! 



5.3 Changing the output ages of the spectra 

ages.dat contains the ages in Myr (one per line, integer) at which the spectra are printed. You may 
change these data (do not forget to type <return> at the end of the file). If necessary, modify the 
parameter nmaxtimesimpr at the beginning of spectra, f (maximal number of printed spectra). 

5.4 Introducing other filters 

You may include other filters (see section 2.3.3): 

• Change the number of filters on the first line of filters.dat. 

• At the end of the file, write on the same line (with blanks between them) : 

— the number of wavelengths defining the passband of the filter; 

— the type of transmission; 

— the type of calibration; 

— the name between quotes; 

— comments (optional). 

• Write each wavelength (A) and the corresponding transmission on one line. Do not forget to type 
<return> at the end of the file. 

• Run calib to obtain the calibrations of the filters in calib.dat. 

5.5 Printing other quantities 

colors may print other quantities: 

• Normalized absolute magnitudes [$, §] mag(j ,i) (magnitude at time time(j) in the filter number 
i) and derived colors [§]. 

• Normalized "monochromatic" luminosities [f, §] f luxf ilter (j , i) (erg.s^.A -1 ) or their ratio to 
the solar luminosity in the filter (f luxf ilter (j ,i)/fluxsol(i)). 

• Normalized luminosity [f, §] Lumline(j ,i) (erg.s -1 ) of the emission line number i (see HII.dat) 
at time time(j) or its equivalent width [§] EW(j ,i) (A). The data for the nebular lines are given 
in HII.dat at lines 84 to 144. Each line contains the wavelength of the emission line, the ratio of 
its intensity to H/3, the name and, finally, the index i used in color s.f. 

To do this, add new lines in colors. f before the instruction close (50) in the following way: 
do j=l,ntimes 

write (50 , format ) variablel (j) , variab le2(j) .. . 
end do 
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